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A TEST FOR RADIAL MIXING USING LOCAL STAR SAMPLES 



(N 

o 

(N 
>Y 



< 

o 

Oh: 
6 



o 
o 
t> 

O 

(N 



13 



JlNCHENG YU 

Key Laboratory for Research in Galaxies and Cosmology, Shanghai Astronomical Observatory, 

Chinese Academy of Sciences, 80 Nandan Road, Shanghai 200030, China 
Graduate University of the Chinese Academy of Sciences, 19A Yuquanlu, Beijing 100049, China 

yujc@shao.ac.cn 

J. A. Sellwood and Carlton Pryor 

Department of Physics and Astronomy, Rutgers University, 
136 Frelinghuysen Road, Piscataway, NJ 08854 
sell wood, pryor@physics.rutgers.edu 

Li Chen and Jinliang Hou 

Key Laboratory for Research in Galaxies and Cosmology, Shanghai Astronomical Observatory, Chinese 
Academy of Sciences, 80 Nandan Road, Shanghai 200030, China 
chenli, houjl@shao.ac.cn 
Accepted to appear in ApJ May 29, 2012, first submitted Sep 30, 2011 

ABSTRACT 

We use samples of local main-sequence stars to show that the radial gradient of [Fe/H] in the thin disk 
of the Milky Way decreases with mean effective stellar temperature. Many of these stars are visiting 
the solar neighborhood from the inner and outer Galaxy. We use the angular momentum of each star 
about the Galactic center to determine the guiding center radius and to eliminate the effects of epicyclic 
motion, which would otherwise blur the estimated gradients. We interpret the effective temperature as a 
proxy for mean age, and conclude that the decreasing gradient is consistent with the predictions of radial 
mixing due to transient spiral patterns. We find some evidence that the trend of decreasing gradient with 
increasing mean age breaks to a constant gradient for samples of stars whose main-sequence life-times 
exceed the likely age of the thin disk. 

Subject headings: Stars: abundances — galaxies: kinematics and dynamics 



1. INTRODUCTION 



A radial abundance gradient in the Milky Way disk to- 
day is w ell-established from analyses of the in terstellar 
medium (|Shaver et ali Il983l : iBalser et aZj|201lh . Metal 
abundances are believed to inc rease over time and most 

Boissier fc Prantzoslll999l: 



chemical evolution models (e 



Hou et ali l2000t IChiappini et a/j|2001l : iNaab fc Ostrikerl 



20061 : Fu et al.ll2009l ) predict that a radial abundance gra- 



dient is created early and may even have been steeper in 
the past. This picture is consistent with "inside-out" disk 
formation, in which the stellar populatio n at increasing 
radii is younger and more metal poor (e.g. Ide Jongj ll999; 
iMunoz-Mateos et aLlfeOOTf ). 

The chemical composition of a star is that of the cloud 
of gas from which it formed. Abundance gradients are 
theref ore also observed in the disk of the M i lky Way in 
stars (INordstrom et alj|2004l ; lHavwoodll2008l : iLuck et ali 
l2011l : [Schlesinger et q/j|2011h. in planetary nebulae C Maciel 



et al. 120071: IStanghcll ini fc Havwoodll20irjf). a nd m open 
star clusters (|F7iel et aLll2002l ; IChen et aLll2003D . While all 
agree that abundances in the thin disk decrease outwards, 
the precise slope differs from sample to sample and element 
measured. 

The orbi ts of disk star s generally become mor e eccentric 
over ti me (IWielenlll977l: iHolmberg et aZ.ll2009l; Aumer & 
Binncy 2009) , which "blurs" measurements of the metal- 
licity gradient when the instantaneous radii of older stars 
are used. Eccentric motion of any amplitude, projected 
to the mid-plane of an axisymmetric potential, can be de- 
scribed as a retrograde epicycle abo ut a guiding center 
that orbits at a constant angular rate ()Binnev fc Tremaind 
2008). The guiding center, or home, radius of a disk star is 



directly related to the z-component of the angular momen- 
tum, L z , about the Galactic center and, since angular mo- 
mentum is conserved, can be determined from any point 
around the epicycle. Thus measurement of the angular 
momenta of stars allows us to eliminate the blurring ef- 
fects of epicyclic motion s and simplifies the inter pretation 
of metallicity gradients. INordstrom et ali (|2004 ) used the 
mean Galacto-centric radius for the same reason, which 
requires integrating the orbit in some adopted potential — 
angular momentum is more direct. 

In the absence of radial mixing, chemical evolution mod- 
els in which the metallicity of the ISM increases at each 
radius over time would predict a tight correlation between 
the metallicity and age for stars of a given home radius. 
This prediction is not supported by the data, and the 
metallicity distributi on of all but the younge st stars shows 
a broad distribution (lEdvardsson et aElll993l: Nordstrom et 
al. 120041 ), which persists even after correcting for epicyclic 
blurring (Ha vrooll2l)08h. 

Sellwood & Binney (|2002f ) argued that stars in a disk 
galaxy are shuffled in radius over time by the effects of 
transient spiral arms. They showed that a star near the 
corotation resonance of a spiral could gain or lose enough 
angular momentum for its home radius to change by up 
to 2 kpc, leaving the star at its new mean galacto-centric 
distance with no increase in its epicyclic amplitude. The 
resulting radial mixing is described as "chur ning" and their 
result h as been confirmed and extended dRoskar et ali 
2008 aHbl. ILoebman et al.ll2011l:lMinchev et aLll2011l: Bird et 



al. 2012 ; ISolwav et al. 2012f ). The broad spread of metal- 
licities among stars having common ages and home radii 
today is then naturally explained as reflecting the past gra- 
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dient of metallicity o f the I SM at their various birth radii. 
Sch onrich fc Binnevl (|2009t ) provide a model for chemical 
evolution of t he Milky Way that wa s the first to include 
radial mixing. iLoebman et all (|2011[ ) find good agreement 
between the predictions of their simulations, that manifest 
substantial radial mixi ng, with the meta ll icity distribution 
of stars from the SDSS (jlvezic et aLll2008f ); lLee et all (|201lf ) 
find other evidence supporting radial migration. 

An additional consequence of radial mixing, or churn- 
ing, is that the radial metallicity gradient of a generation 
of stars is gradually flattened over time. Nordstrom et 
al. (|2004f ) and iHavwoodl (|2008l ), from local stars suggest 
that the gradient may become shallower with increasing 
age. ICasagrande et all (|201lD report little evolution of 
the radial metallicity gradient over the past ~ 5 Gyr, us- 
ing infrared magnitudes and improved models to revise 
the ag e estimates for the stars in the iNordstrom et al\ 
( 2004 ) sample. Estim ating the age of an individual star 
(see ISoderbloml |2010| . for a recent review) is a delicate 
art and the results even for well-observed mai n-sequence 
stars can be contentious (e.g. iReid et aZ.1 120071; Holmberg 
et al I2007T1. F urthermore, the eviden ce from star clus- 
ters (|Friel et aLll2002l ; IChen et a/.ll2003D suggests that the 
radial metallicity gradient may steepen wit h age, in line 
with i nside-out models of galaxy formation. iMaciel et all 
(2001]), from a study of planetary nebulae abundance data, 
also find evidence for a steeper gradient amo n g thei r older 
objects, although IStanghellini fe Havwoodl (|2010f ) reach 
the opposite conclusion that older PNe show a shallower 
gradient. We comment on these apparently conflicting re- 
sults in §4.2. 

Here we study the age-dependence of the metallicity gra- 
dient using a sample of stars in the solar neighborhood. 
Instead of estimating ages of individual stars, we adopt 
the effective temperature of a collection of main sequence 
stars as a proxy for their mean age, in the same man- 
ner that lAumer fc Binnevl (|2009t ) used color. Since the 
main-sequence lifetime of a star is shorter for hotter stars, 
the mean ages of main sequence disk stars grouped by ef- 
fective temperature must be lower for the hotter groups. 
(Of course, this argument requires the reasonable assump- 
tions that stars have been forming with an approximately 
constant IMF and at a roughly uniform rate over the life- 
time of the Milky Way disk.) The trend of increasing 
mean age with decreasing temperature must change to 
a constant at the point at which the main-sequence life- 
times of stars exceed the age of the disk. Note that the 
main-sequence turn-off for stars of age 10 Gyr occurs for 
5 500 < T eff < 6 000 K over the range 0.4 > [Fe/H] > -0.5 
(|Demarque et aZ.ll2004l V 

As explained above, we eliminate epicyclic blurring by 
estimating the specific angular momentum of each star, 
which requires knowledge of its full 6D phase space coor- 
dinates: i.e. distance and proper motion, as well as radial 
velocity and sky position. A radial velocity can be mea- 
sured spectroscopically at any distance, but a reliable dis- 
tance can be obtained only for nearby stars. Furthermore, 
distance uncertainties factor into velocity components in 
the sky plane adding to the desirability of restricting at- 
tention to stars close to the Sun. Also, by focusing on 
nearby stars, any variations in L z across our sampling vol- 
ume due to possible departures from axisymmetry in the 



Galactic potential will be negligible. 

A benefit of stellar epicycle excursions is that they bring 
stars to the solar neighborhood. Thus a sample of very 
local stars will span a significant range of specific angular 
momenta and, therefore, home radii. 

In this paper, we assemble a sample of local, main- 
sequence stars having estimated effective temperatures, 
metallicities and full space motions. We draw these stars 
from three sources as described in the next section. 

2. STAR SAMPLES 
2.1. Geneva-Copenhagen sample 

The Gene va-Copenhagen survey of nearby stars ( Nord- 
strom et al. I2004J ) used Hivvarcos positions and proper 
motions, supplemented by spectroscopic observations, to 
construct a homogeneous sample of nearby mostly F and 
G dwarf s tars. W e employ their updated catalog (Holm- 
berg et al. 120091 ) that uses their revised temperature cali- 
bration and the impr oved astrometr y from the reanalysis 
of Hipparcos data by Ivan Leeuwenl (|2007f ). This monu- 
mental effort has resulted in a large sample of stars in 
the sol ar neighborhood with d istances and full space mo- 
tions. iHolmberg et al\ ((2009) do not supply individual 
uncertainties for each velocity, but assert that they are 
believed accurate to 1.5 km s -1 , with the greatest contri- 
bution coming from distance uncertainties. 

The machine readabl e table of 16 682 stars made avail- 
able bv lHolmberg et al\ (|2009f ) includes sky positions, dis- 
tances, and the (U, V, W) components of the star's motion 
relative to the Sun in Galactic coordinates^ We discard 
those having no distance, no (U, V, W) velocities, or no 
estimate of T e g. We do not use the disputed ag e esti- 
mates fe.q. lReid erZll2007l ; lHolmberg e£ kl 120071 ) in the 
present paper. We correct for Solar motion relative to 
the local standard of rest ( LSR) by adding (U, V , W)q = 
(11.10,12.24,7.25) km s" 1 (jSchonrich et allhoiCt ) to the 
tabulated velocitie s. 

Nordstrom et al. (|2004D note that they obtained [Fe/H] > 
0.4 for a small number of stars, though none is within 40 pc 
of the Sun. Furthermore, these high values are found only 
for stars with T e g > 6200 K. They argue that the coinci- 
dence of high T c ff and [Fe/H], together with other informa- 
tion, suggests that the extinction to these stars has been 
overestimated. We therefore discard a further 29 stars 
having [Fe/H] > 0.4. 

In order to select nearby disk stars, we further restrict 
the sample to stars whose best estimate of the distance 
is within 500 pc, \U\, \V\ < 80 km s _1 , and retain only 
those having an energy of vertical motion about the Galac- 
tic mid-plane, E z = 0.5(z 2 i/ 2 + W 2 ) < 392 (km s" 1 ) 2 , with 
the verti cal frequency v = 0.07 km s -1 pc -1 ( Binnev & 
Trcmaine l2008l ). giving them a maximum vertical excur- 
sion of ±400 pc. Thus we select only stars that have a high 
probability o f being thin disk stars, in the same spirit as 
Bcn sbv et all (|2003l ). but not in exactly the same manneiQ 
(see §3]for further discussion of selection effects). Our final 

1 These Cartesian velocity components are oriented such that U is 
towards the Galactic centre, V is in the direction of Galactic rotation, 
and W is towards the north Galactic pole. 

2 Adopting their selection criterion for thin disk stars, changes the 
number of selected stars by a few hundred, but our conclusions are 
unaffected by this marginal revision. 
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Fig. 1. — The densities of stars in the space of T e fj and [Fe/H] in 
our selected GCS (red) and RAVE (blue) subsamples. 

sample contains 11 877 stars that we use in this analysis. 

The median distance of the selected stars is 75 pc and 
the sample within 40 pc of the Sun is believed to be near 
complete. We adopt Ro — 8 kpc, the circular speed of the 
LSR to be Vo = 220 km s _1 , and compute the Galacto- 
centric angular momentum, L z , using the in-plane distance 
of the star from the Sun and the peculiar velocities (U, V) 
corrected for the Sun's motion. 

The red contours in Fig. Q] show the distribution of the 
selected GCS stars in the space of T g and [Fe/H]. The 
range of T e ff indicates that the sample extends outside 
the spectral typ es F and G, reflec t ing th e selection crite- 
ria described by iNordstrom et all (|2004f ). These authors 
estimate uncertainti e s of < 0.1 in [Fe/H] and 94 K in T e g . 

Casagrande et al. (|2011[ ) present a recalibration of T e g 
and [Fe/H] for the GCS stars that folds in the infrared flux 
from a star. They obtain slightl y higher values for T e g 
and [Fe/H] than those derived bv lHohnberg et al] |2009) 
and they too attempt to assign ages for the individual 
st ars. The ad ditional informat ion used by Casagrande et 
al. (|2011l ) should lead to improved values for the derived 
quantities, but they claim reliable results for fewer than 
half the GCS stars (their irfm sample) and for none of the 
cooler stars. We therefore co ntinue to use the calibr ation 
for the whole GCS sample bv lHolmberg et all (|2009l ). 

2.2. RAVE sample 

The RAVE survey (jSteinmetz et aLll2006t ) plans to mea- 
sure the heliocentric radial velocities and stellar parame- 
ters for about a million stars in the southern sky having 
apparent magnitudes in the range 9 < mi < 12; the first 
77 461 are available in the third da ta release (Siebert et 
al. 120111 ). The typical uncertainty in the radial velocity 
is < 2 km s _1 , but the distance to most stars has to be 
judged photometrically and most proper motions are from 
ground-based data. Thus three of the phase space coordi- 
nates for each star are of much lower quality than are those 
in the GCS, although this weakness will be compensated 
by a much larger sample size when the survey is complete. 

We have downloaded the on-line table of the third data 
release from the RAVE website and selected a subset of 
stars for analysis^ We estimate distances to these stars 

3 The paper by Coskuno glu et al\ 1120111) presenting an analysis of 



by fit ting to the Yonsei-Yale isochrones (jDemarque et all 
using a method rela ted to that described by Breddels 
et al. (|2010l . see also Burnett et al. 2011 for an alternative 
method). We adopt many of their selection criteria: we 
require the spectral signal-to-noise parameter S2N > 20 
with a blank spectral warning flag field; the parameters 
[M/H], log(g), and T c g to be determined; and the stars to 
have J and K s magnitudes from 2MASS with no warning 
flags about the identification of the star or the 2MASS 
photometry Unlike those authors, however, we have kept 
stars with b < 25° on the grounds that extinction for the 
nearby stars that interest us will not be large enough in the 
near IR to severely bias our distance estimates. As we are 
here interested only in nearby main-sequence stars that are 
members of the disk population, we also make preliminary 
cuts to eliminate stars with log (g) < 4, T c g > 10 4 K, and 
with |iv| > 120 km s . We have treated the tabulated 
[M/H] a s equivalent to [Fe/H1 since the corrections ( Zwit- 
ter et al. 200§D are generally within the uncertainties. We 
also show in §2.3 that the [M/H] values for RAVE stars 
agree with the [Fe/H] values for the stars in common with 
the GCS. 

We estimate the absolute J magnitude of each selected 
star by matching the estimated [Fe/H], log(g), T e g, and 
J-K s color to values in the isochrone tables for stars of all 
ages and all values of [a/Fe], rejecting a few more stars 
for which the best match % 2 > 6. We consider the closest 
match in the tables to the given input parameters to yield 
the best estimate of the absolute magnitude from which 
we estimate a photometric distance using the apparent J- 
band magnitude|f| We save the values of [Fe/H], log(g), 
T e ff, and J-K s color of the closest matching model star in 
the isochrone table; Monte-Carlo variation of the stellar 
parameters about this saved set of values suggests that 
distances have a fractional precision of 30% - 50%, with 
some larger uncertainties. 

We use the proper motions in equat orial coordinate s 
tabulated in RAVE, mostly from Tycho-2 (|H0g et aZ.ll2000h FI 
which we then combine with the radial velocity and po- 
sition to determine the heliocentric velocity in Gala ctic 
componen ts (U.V.W) dJohnson k. Soderbloml Il987l; Pi- 
atek et al. I2002T ). We estimate uncertainties in these ve- 
locities from 500 Monte Carlo re-selections of all the stel- 
lar parameters that affect the distance estimate, adopt- 
ing cr(J) = 0.03 mag, cr(J - K s ) = 0.042 mag, cr(T off ) = 
300 K, g(log q) = .3 dex, and <r([Fe/H]) = 0.25 dex 
(jBreddels et aLlr2010D . as well as the tabulated radial ve- 
locity and proper motion uncertainties. 

We exclude stars more distant than 500 pc, correct the 
space velocities for solar motion and apply the same re- 
strictions as for the GCS sample to select only those with 
a high probability of being thin disk stars. 

Fig.[5]shows histograms of distances and of uncertainties 
in L z for the 9,803 remaining stars. The median distance 

the radial abundance gradient from these same stars appeared while 
this paper was in review. 

4 Zwittcr et al. (2010) describe a similar method, but define a "most 
likely" estimate of the absolute magnitude that differs from our 
"best" estimate. The difference is likely to be well within the uncer- 
tainties for the main-sequence stars considered here. 

5 We have not used the newly available proper moti ons from the 
fourth release of the Southern Proper Motion Catalog llGirard et al\ 
I20TID 



4 



Yu et al. 



400 



200 




200 400 

dist (pc) 




2 4 6 8 10 

(%) 



400 



200 




200 400 

dist (pc) 



2 4 6 8 10 

<KO/l (%) 



Fig. 2. — Histograms of distances and L z uncertainties for the 
selected sample of 9803 RAVE stars. 



Fig. 3. — Histograms of distances and L z uncertainties for the 
selected sample of 11 022 M-dwarf stars. 



is 220 pc from the Sun and uncertainties in L z are generally 
smaller than 5%. 

The blue contours in Fig. [T] show the distribution of 
selected RAVE stars. The distribution in [Fe/H] values is 
slightly broader than for the GCS stars, reflecting in part 
at least the lower precision of the spectra. However, the 
rms scatter in [Fe/H] appears to be ~ 0.2, sugg e sting the 
uncertainty of 0.25 estimated by iSiebert et ali (|2011h is 
somewhat pessimistic for this sub-sample of main-sequence 
stars. These authors also suggest uncertainties of 300 K in 
their estimates of T e ff, and the broader range of T e g than 
in the GCS sample is therefore real. 

2.3. Stars in common 

The overlap between the GCS and RAVE samples of stars 
is quite small but, of those we have selected, precisely 
100 stars are positionally coincident. This small degree 
of overlap allows us to compare the spectral parameters, 
T e ff, [Fe/H], and radial velocity (v r ), derived separately by 
the two teams from their independent spectra. 

We find that differences in the estimated v r for the same 
stars average just 0.32 km s _1 higher in RAVE than in 
GCS, with an rms scatter about this mean of only 2.3 km s _1 . 
Similarly, differences in T e g have a mean of 203 K, again 
higher in RAVE than in GCS, with an rms scatter about 
the mean of 202 K. Finally, the mean difference in [Fc/H] 
is 0.057, lower in RAVE than in GCS, with an rms scatter 
of 0.18. We have looked in vain for systematic variations 
of these differences, which appear to be consistent with 
random scatter. 

Thus all three spectral parameters agree between the 
independent measures to an accuracy that is significantly 
better than expected from the quoted uncertainties, indi- 
cating that both teams have been very careful with their 
calibrations and conservative with their uncertainty esti- 
mates. 

2.4. M-dwarf sample 

Th e SDSS (jYork et a/.ll2000f ) and Segue2 (jYannv et all 
2009]) surveys are complete, but sample fainter stars than 
RAVE (the magnitude range for Segue2 was 14.0 < g < 
20.3) that are therefore generally m ore distant. Sin ce the 
Sloan spectral parameters pipeline ( Lee et alll2008h does 



not attempt to fit stars with T e g < 4500 K, almost all 
main sequence stars with estimated parameters in SDSS 
DR8 are more distant than 500 pc . 

Fortunately, IWest et~al\ (|201lD provided a catalog of 
70 841 M-dwarfs from SDSS DR7. Their table provides a 
photometric estimate of the distance to each star, as well 
as a spectroscopic radial velocity, and proper moti on from 
USNO-B/SDSS catalog dMunn et aZJ l2004L l2008f>. They 
suggest that distance uncertainties are typically about 20% 
and uncertainties in radial velocity are 7-10 km s _1 . None 
of these stars are in RAVE (because they are in different 
parts of the sky) and all are much cooler than the GCS 
stars. 

As lWest et ali (| 2 1 ll ) recommend, we selected stars with 
the "goodPM" and "goodPhot" flags set to 'true', and the 
"WDM" flag set to 'false' to eliminate possible binaries 
with a white dwarf companion, which reduces the sample 
to 39 151 stars. We further excluded stars having no ra- 
dial velocity as well as those with no distance estimate or 
for which the estimated distance exceeded 500 pc. As for 
the GCS and RAVE stars, we correct the space velocities 
for solar motion and apply the same restrictions to select 
only those with a high probability of being thin disk stars, 
leaving us with a final sample of 11 022 stars. 

We estimated uncertainties in the velocity components, 
(U, V, W) by combining the 20% distance uncertainty, a 
10 km s _1 uncertainty in the radial velocity, and the proper 
motion uncertainties. The distribution of distances and 
fractional L z uncertainties is shown in Fig. [3] many of 
these intrinsically faint stars lie within 200 pc and, again, 
angular momentum uncertainties are typically < 5%. 

As estimates of T c g or [Fe/H] are not easily derived from 
low resolution spectra of such cool stars. IWest et ali (|201ll ) 
do not attempt to provide these quantities in their cata- 
log. Instead, they provide spectral sub-class, which can be 
related to T e ff, and a titanium oxide index, CnOi which is 
believed to be an indicator of metallicity. 

3. RESULTS 

3.1. GCS and RAVE stars 

Figs. show the metallicity vs. angular momentum 

of GCS and RAVE stars each separately binned into various 
temperature ranges. The bins are chosen so as to have 
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Fig. 4. — The distributions of GCS stars in the space of L z and 
[Fe/H] when grouped by effective temperature. The temperature 
bins were chosen to contain equal numbers of stars and the average 
temperature of the stars in each bin is given. The large open circles, 
which are purely to illustrate the trends, show the median [Fe/H] 
for stars with 1300 < L z < 2100 divided into six bins in each panel. 
The unit of angular momentum is km s — 1 kpc. 

equal numbers of stars, and the average T c g of the stars in 
each panel is shown. The density of points in each panel 
is indicated by the contours. 

If the rotation curve of the Milky Way were flat at a 
constant circular speed of 220 km s _1 , then stars with 
L z — 1 210 km s _1 kpc would have home radii of 5.5 kpc, 
while the home radii of those with L z = 2 200 km s _1 kpc 
would be 10 kpc. (Note that these numerical values will 
simply scale with our choices of Rq and Vq.) Thus, al- 
though the stars in our sample are passing close to the 
Sun right now, many are visiting from quite far afield. 

A number of trends can be seen: The spread in L z is 
greater for the cooler stars, reflecting the fact that the 
higher velocity dispersions of older stars allows stars hav- 
ing a wider range of L z to visit the solar neighborhood. 
The rms scatter in [Fe/H] is in the range 0.19 < a < 
0.23 for the RAVE stars in each panel, which is less than 
the nominal uncertainty, as noted in W2.2\ The spread 
in [Fe/H] for the GCS stars is generally smaller, rising 
from 0.14 in the second panel to 0.21 for the coolest stars. 
The spread for the hottest GCS stars (0.16) is anomalously 
high, due to incomplete removal of stars that were over- 
corrected for reddening, as noted in ij2.ll While the large 
spread in [Fe/H] is partly due to uncertainties and partly 
real, each panel has a weak declining trend of metallicity 
with increasing L z that reflects the metallicity gradient in 
the Galaxy. We have added the open circles, which show 
the variation of the median in bins of L z in each panel, 
simply to make this last trend more clearly visible. 

3.2. Trend in gradient 




-0.8 1 1 1 1 1 1 1 ' — 1 1 1 1 1 1 ' — ' 

1200 1600 2000 2400 1200 1600 2000 2400 

L z (km/s kpc) L z (km/s kpc) 



Fig. 5. — Same as Fig. [4] but for the RAVE sample. 



We adopt Theil's non-p arametric method (Hollander & 
Wolfe Il999f ) to quantify the metallicity gradient in each 
panel, which we plot as filled symbols in Fig. [SJ with the 
dashed error bars showing the 95% confidence limits. This 
method returns the median slope between every pair of 
points in the sample with the confidence limits being de- 
termined by the range of the pairwise values. We have con- 
verted the gradient with angular momentum to the more 
usual dex per kpc by multiplying by our adopted circular 
speed Vb = 220 km s -1 . Our abundance gradient esti- 
mates therefore assume only a locally flat rotation curve 
and a distance scale set by adopting Rq — 8 kpc. 

While these estimates of the metallicity gradient are de- 
rived from stars that happen to be close to the Sun at the 
present time, the sample includes stars from a broad range 
of home radii, as noted above. Stars whose ages greatly 
exceed the radial epicycle period, ~ 170 Myr at the so- 
lar radius, should be uniformly distributed in radial phase 
around their epicycles. Thus both the spread in metallici- 
ties of the stars in our sample and the gradient we measure 
are characteristic of the stars at their home radii, and our 
estimates of the metallicity gradient are based on a radial 
range of several kiloparsec. 

The apparent flattening of the metallicity gradient with 
decreasing T c g is suggested by the RAVE stars (left panel 
of Fig. [5]) and is more convincing in the better-quality 
data of the GCS sample (right panel). As this is our main 
result, we use Monte-Carlo simulations to check whether 
the uncertainties in the slopes are small enough to justify 
our claim of a trend. We generate 100 new samples, by 
randomly resampling stars, adding normally distributed 
values to T e ff, [Fe/H], and L z to simulate measurement er- 
rors. The size of the errors are determined by the appro- 
priate uncertainties from each survey as described above 
or, for the L z values of the RAVE stars, by propagating 
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Fig. 7. — TiO index, C/riO vs - angular momentum for the M- 
dwarf sample, divided by spectral sub-type into two roughly equally 
numerous groups. 



Fig. 6. — The metallicity gradient, d[Fe/H]/dL z estimated from 
the GCS sample shown in Fig. [4] (right panel), and from the RAVE 
sample in Fig. \E\ (left panel). The filled symbols show the best esti- 
mate from Theil's algorithm, with 95% confidence limits (dashed), 
the crosses with solid error bars show the standard errors derived 
from a Monte Carlo simulation with the adopted uncertainties. 



our own estimates of the distance, proper motion and ra- 
dial velocity uncertainties. We divide each sample into 
equal numbers of stars in successive bins of T G ff and com- 
pute the slopes. The means and standard errors of these 
re-estimates are also shown by the crosses and solid er- 
ror bars in Fig. [6] Note that the solid error bars from 
our Monte-Carlo simulations indicate the 1-a dispersion, 
whereas the dashed error bars from the Theil algorithm 
show the 95% confidence interval; the two independent es- 
timates are therefore consistent, and the flattening of the 
slope with decreasing T e g seems to be confirmed. Note 
also that the somewhat shallower trend in the Monte- 
Carlo re-estimates from the RAVE data is expected, since 
resampling the data in this manner inserts additional un- 
certainties over and above those already present in the 
data, which inevitably weakens a trend. This effect is 
more marked in the RAVE sample because uncertainties 
in their data are larger. 

It also occurred to us that the change of slope could 
perhaps be a systematic effect due to the above-noted in- 
creasing spread in L z with decreasing T c g - i.e., the greater 
spread in L z could simply make the fitted slope appear 
shallower, while in fact the data may be consistent with a 
constant slope. We tested for this possibility by creating 
some pseudo star samples that adopt a linear relation be- 
tween [Fe/H] and L z that is the same for all the stars in 
either the RAVE or the GCS sample, irrespective of their 
temperature, which is the null hypothesis that we argue 
Fig. |5] disproves. For each star in the sample, we use the 
measured L z to assign a new [Fe/H] from a Gaussian dis- 
tribution about the mean from the adopted trend with a 
dispersion equal to that about the linear regression line for 
the entire sample. We then bin this pseudo-sample by T e g 
as before and test whether the cooler stars, that have the 
larger spread in L z , have a shallower apparent slope. With 
500 realizations for each of the GCS and RAVE samples, 
we confirmed that we recover the original input slope in 
every bin. The absence of any systematic trend with T e g 
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Fig. 8. — The estimated metallicity gradient d[Fc/H]/dL z from 
the GCS stars (red) and from the RAVE stars (blue) reproduced 
from Fig. [6] The green points show d^io/dLz from the M-dwarf 
stars in Fig. [7] further subdivided by spectral sub-class. The error 
bars show 95% confidence limits from Theil's algorithm. 



in this test increases our confidence that the trend we find 
in the data in Fig. [5] is real. 

Furthermore, both samples also hint at a break in the 
trend with temperature around T e g < 6 000 K, although 
these points alone are also consistent with a uniform trend. 
In order to find extra evidence for this possible break of 
slope, we attempt to include the M-dwarf stars into our 
analysis. 

3.3. M-dwarf sample 

Fig. [7] shows the distribution of L z and (riO for the M- 
dwarf stars separated into two groups by spectral class. 
This metallicity indicator again has a noticeable gradient 
with L z . 

As for the hotter stars, we have used Theil's algorithm 



Radial mixing 



7 



to estimate the gradient dCTio/dL Zl dividing the M-dwarf 
sample into four groups by stellar sub-class, and plot the 
results in green in Fig. [5] (The abscissae were determined 
from the relation T c g = 3700— 1505" K, where S is the spec- 
tr al sub-class, fitted by eve to the compilation by Revle et 
al. (|201lf ).') The M-dwarfs have a constant metallicitv gra- 
dient that appears independent of T c s . 

Because they are measured from a different quantity, 
the values of dQno/dL z cannot be compared directly with 
those of d[Fe/H]/dL z from the hotter stars, shown in red 
and blue. We tried using the e mpirical relat i on be tween 
Ctio and [Fe/H] suggested by IWoolf et all (|2009l ), but 
many [Fe/H] values required an extrapolation beyond the 
range calibrated by their study, and the resulting mean 
[Fc/H] is unreasonably high. Thus we lack a direct basis 
to relate the green points in Fig. [3] to the others. 

However, we note that if £riO is proportional to [Fc/H], 
the gradient d[Fc/H]/dL 2 will scale with dQrio/dL Zl and 
we would merely need to shift the green points vertically 
by some fixed factor to bring them onto the same scale. 
Since an anti-correlation between CtiO and [Fe/H] seems 
unlikely, the green points probably could not be shifted 
above d[Fe/R]/dL z = 0. Thus the data from the M- 
dwarf stars do seem to confirm that the decreasing trend 
of metallicity gradient with T Q g among the hotter stars re- 
ally does break to no trend around T e g < 6 000 K, i.e. for 
stars whose main-sequence lifetimes arc ~ 10 10 yr. 

4. DISCUSSION 



Ma ny authors (e.g. Lee et all 1201 ll ; iCasagrande et aj\ 



120111 : ICoskunoglu et aZ.ll2011h express concerns that est 



mates of the metallicity gradient from a sample of stars de- 
pend upon the adopted selection criteria. This is because 
metallicity is a function of age and, relative to younger 
stars, older stars typically have larger peculiar velocities, 
which causes an asymmetric drift and greater vertical mo- 
tion. This well-known kinematic trend must cause an age 
bias in kinematically selected samples. 

4.1. Possible selection effects 

In this work, we apply kinematic and vertical motion 
cuts to eliminate most thick disk and halo stars, which may 
well also eliminate some old thin disk stars that should 
perhaps be included in our samples. However, any bias 
in our analysis created by excluding older stars because of 
their higher peculiar velocities will cause lower tempera- 
ture bins, which are the only bins to include stars of the 
old disk, to have unrepresentatively few older disk stars. 

Radial mixing implies that these possibly-missing stars 
would be members of the most well-mixed population with 
the lowest radial metallicity gradient, and adding signifi- 
cant numbers of them to the lower T c g bins would only 
increase the trend we observe. In fact, we find that the 
trend of metallicity gradient with T e g is little affected by 
limiting the vertical excursions of selected stars to either 
300 pc or 500 pc, suggesting that this cut is not biasing 
our result. Further restricting the vertical oscillations to 
< 300 pc, on the other hand, does exclude a dispropor- 
tionate fraction of older stars and the trend with T e ff is 
weakened, as our argument predicts. Naturally also, the 
number of selected stars decreases with a more severe re- 
striction increasing the statistical uncertainty in our esti- 



mate of the metallicity gradient. 

We do not use any measure of chemical composition to 
select stars, and our sample may therefore include some 
stars with enhanced [cc/Fe]. Whether such abundances 
best ch aracterize thick disk stars is the sub ject of current 
debate (jlvezic et al\ 120081 : iLee et aZj|201ll ) while the ex- 
istence of a thick disk that is disti nct from the thin has 
been questioned (|Bovv et q.Z.1 1201 H ) . Our study is unaf- 
fected by these questions, because we are interested only 
in how the metallicity gradient changes with the mean age 
of the stars and, provided we have excluded halo stars 
that are little affected by spiral activity, we do not need 
to consider whether each star m ay or may not belong to a 
particular population. Note that ISolwav et al\ (|2012f ) have 
shown that mixing is only gradually weakened by increased 
random motion and vertical thickness. 

A separate issue is that all the stars in this study are 
close to the Sun at the present time, and we therefore 
cannot include stars with home radii far from the Sun but 
with radial excursions that arc too small to bring them into 
the solar neighborhood. However, provided the stars we 
do see have abundances that are representative of those of 
their T Q g and home radius, the omitted stars will not affect 
our estimate of the slope. Since we are not selecting on 
metallicity, there should be no bias in our slope estimates. 

4.2. Comparison with other work 

Casagrande et al. (|201ll ) re-analysed the GCS using a 
revised calibration of T e g and [Fe/H] and new estimates of 
the individual stellar ages. We have tried to compare our 
results with theirs, which is not straightforward because 
our bins of T e ff each contain stars havi ng a broad range 
of ag es. We have used isochrone tables (|Demarque et al\ 
2004) to estimate the main-sequence lifetimes of the indi- 
vidual stars in the GCS sample, using the T c g and [Fe/H] 
of each, and so find the average main-sequence lifetime of 
the stars in our separate T Q g bins. The mean maximum 
age of stars defined in this way increases from about 2 Gyr 
for our hottest bin to 10 Gyr for the coolest. Using these 
mean maximum ages to replot the trend of metallicity gra- 
dient in our Fig. 8, we find reasonable agreement with the 
gradient valu es as a cumulativ e funct ion of age shown in 
Fig. 18(d) of lCasagrande et al\ (|2011h . 

Our estimated gradient is in t he range of other m e asure- 
ments from stars. For example, ICoskunoglu et ali (| 2 1 lh 
derived a metallicity gradient of -0.051 dex kpc -1 for F- 
dwarfs, and they also find a shallower gradient for the 
slightly older G-dwarfs. While their sample clearly over- 
laps with ours, they adopted different selection criteria 
and employed a different technique to determine the home 
radii, it is encouraging that they reached similar conclu- 
sions ^UsingcUsk^Cerjhejd^whic^ 



et al. (J2011J) derived [Fe/H]= -0.055i? + 0.475, where R 
is the Galacto-centric radius, again in tolerable agreement 
with our grad i ent es timate for the hotter stars. 

Chen et al. (|2003T l estimate a radial gradient of [Fc/H] 
of -0.063 ± 0.008 dex kpc" 1 from 118 star clusters. Both 
they and iFriel et "all (|2002l ) report that, when the clus- 
ters are divided by age, the slope appears steeper for the 
older sample, which is the opposite of our finding here. 
Note that the radial distribution of star clusters in these 
papers extends outwards from the solar radius, whereas 



8 



Yu et al. 



our sample of stars overlaps the solar radius and extends 
inwards to ~ 5 kpc - thus the different variations with 
age could indicate the absence of mixing in the outer disk. 
However, other systematic differences may complicate this 
interpretation; e.g., the outer di s k ope n clusters in the 
sample assembled bv lChen et al\ (|2003[ ) are farther from 
the mid-plane t han a re those near the solar radius. 

Maciel et al. (|2006l ) also claim a steeper abundance gra- 
dient for older planetary nebulae, although their estimated 
gradients for the different elements show a lot of scatter 
and ages are less reliable than for star cluste rs. It should 
be noted that IStanghellini &: Havwoodl (|2010D reached the 
opposite conclusion finding, for a sample of PNe concen- 
trated over the radial range 3 < R < 12 kpc, a shallower 
abundance gradient for the older PNe, which is in better 
agreement with our results. 

5. CONCLUSIONS 

We have used samples of local thin-disk, main-sequence 
stars grouped by effective temperature to demonstrate that 
the radial abundance gradient in the disk of the Milky Way 
flattens as the T c g of the stars in the sample decreases. We 
eliminate the blurring effect of cpicyclic motion by esti- 
mating the angular momentum, L z , of each star about the 
center of the Galaxy, since the home, or guiding center, 
radius of a star is determined only by its angular momen- 
tum. Although the stars in our sample are all within 500 
pc of the Sun, their home radii are spread over a wide 
swath of the Milky Way disk. 

We interpret T e s as a proxy for mean age, and con- 
clude that a shallower gradient for stars of a greater mean 
age. We find some evidence that the trend of decreasing 
metallicity gradient with increasing age breaks to a fiat 
trend around 5 500 < T c g < 6 000 K, for which the main- 
sequence lifetime of a star is approximately 10 Gyr, or 
roughly the age of the thin disk. Flattening of the metal- 
licity gradient is consistent with the predictions of radial 
mixing, or churning, due to spiral patterns in the disk. 

We thank Birgitta Nordstrom for helpful correspondence 
about the apparently metal rich stars in the GCS and an 
anonymous referee for comments that have helped us to 
improve the paper. This work was supported by NSF 
grant AST 09-37523 (PI Newberg), and grant 11061120454 
(PI Deng) from the National Science Foundation of China 
(NSFC), the NSFC grants 11073038 (PI Chen), 11173044 
(PI Hou), Key Project 10833005 (PI Hou), and by the 
Group Innovation Project No. 10821302. 



REFERENCES 

Aumer, M. & Binney, J. J. 2009, MNRAS, 397, 1286 
Balscr, D. S., Rood, R. T., Bania. T. M. & Anderson, L. D. 2011, 
ApJ, 738, 27 

Bensby, T., Feltzing, S. & Lundstrom, I. 2003, A&A, 410, 527 
Binney, J. & Tremaine, S. 2008, Galactic Dynamics 2nd Ed. 

(Princeton: Princeton University Press) 
Bird, J. C, Kazantzidis, S. & Weinberg, D. H. 2012, MNRAS, 420, 

913 

Breddels, M. A., Smith, M. C, Helmi, A., et al. 2010, A&A, 511, 
A90 

Boissier, S. & Prantzos, N. 1999, MNRAS, 307, 857 
Bovy, J., Rix, H-W. & Hogg, D. W. 2011, arXiv: 1111.6585 



Casagrandc, L., Schonrich, R., Asplund, M., et al. 2011, A&A, 530, 
A138 

Chen, L., Hou, J. L. & Wang, J. J. 2003, AJ, 125, 1397 
Chiappini, C, Matteucci, F. & Romano, D. 2001, ApJ, 554, 1044 
Co§kunoglu, B., Ak, S., Bilir, S., et al. 2011, MNRAS, in press 
de Jong, R. S. 1996, A&A, 313, 377 

Demarque, P., Woo, J-H., Kim, Y-C. & Yi, S. K. 2004, ApJS, 155, 
667 

Edvardsson, B., Andersen, B., Gustafsson, B., et al. 1993, A&A, 
275, 101 

Friel, E. D., Janes, K. A., Tavarez, M., et al. 2002, AJ, 124, 2693 
Fu, J., Hou, J. L., Yin, J. & Chang, R. X. 2009, ApJ, 696, 668 
Haywood, M. 2008, MNRAS, 388, 1175 

Hou, J. L., Prantzos, N. & Boissier, S. 2000, A&A, 362, 921 
Girard, T. M., van Altena, W. F., Zacharias, N., et al. 2011, AJ, 
142, 15 

H0g, E., Fabricius, C, Makarov, V. V., et al. 2000, A&A, 355, L27 
Hollander, M. & Wolfe, D. A. 1999, Nonparametric Statistical 

Methods 2nd Ed. (New York: John Wiley) pp. 416-423 
Holmberg, J., Nordstrom, B. & Andersen, J. 2007, A&A, 475. 519 
Holmberg, J., Nordstrom, B. & Andersen, J. 2009, A&A, 501, 941 
Ivezic, Z., Sesar, B., Juric, M., et al. 2008, ApJ, 684, 287 
Johnson, D. R. H. & Soderblom, D. R. 1987, AJ, 93, 864 
Lee, Y. S., Beers, T. C, Sivarani, T., et al. 2008, AJ, 136, 2022 
Lee, Y. S., Beers, T. C, An, D., et al. 2011, ApJ, 738, 187 
Loebman, S. R., Roskar, R., Debattista, V. P., et al. 2011, ApJ. 737, 

8 

Luck, R. E.. Andrievsky, S. M., Kovtyukh, V. V., et al. 2011, AJ, 
142, 51 

Maciel, W. J., Lago, L. G. & Costa, R. D. D. 2006, A&A, 453, 587 
Maciel, W. J., Quireza, C. & Costa, R. D. D. 2007, A&A, 463, L13 
Minchev, I., Famaey, B., Combes, F., et al. 2011, A&A, 527A, 147 
Munoz-Mateos, J. C, Gil de Paz, A., Boissier, S., 2007, ApJ, 658, 
1006 

Munn, J. A., Monet, D. G., Levine, S. E., et al. 2004, AJ, 127, 3034 
Munn, J. A., Monet, D. G., Levine, S. E., et al. 2008, AJ, 136, 895 
Naab, T. & Ostriker, J. P. 2006, MNRAS, 366, 899 
Nordstrom, B., Mayor, M., Andersen, J., et al. 2004, A&A, 418, 989 
Piatek, S., Pryor, C, Olszewski, E. W., et al. 2002, AJ, 124, 3198 
Reid, I. N., Turner, E. L., Turnbull, M. C, et al. 2007, ApJ, 665, 
767 

Reyle, C, Rajpurohit,, A. S., Schultheis, M. & Allard, F. 2011, 
[arXiv:1102. 12631 

Roskar, R., Debattista, V. P., Quinn, T. R., et al. 2008a, ApJL. 684, 
L79 

Roskar, R., Debattista, V. P.. Stinson, G. S., et al. 2008b, ApJL. 
675, L65 

Schlesinger, K. J., Johnson, J. A.. Rockosi, C. M., et al. 2011, ApJ, 

submitted farXiv: 1112.2213} 
Schonrich, R. & Binney, J. 2009, MNRAS, 396, 203 
Schonrich, R., Binney, J. & Dehnen, W. 2010, MNRAS. 403, 829 
Sellwood, J. A. & Binney, J. J. 2002, MNRAS, 336 , 785 
Shaver, P. A., McGee, R. X., Newton, L. M., 1983, MNRAS, 204, 

53 

Siebert, A., Williams, M. E. K., Siviero, A., et al. 2011, AJ, 141, 
187 

Soderblom, D. R. 2010, ARAA, 48, 581 

Solway, M., Sellwood, J. A. & Schonrich, R. 2010, MNRAS, to appear 

(arXiv:1202.1418) 
Stanghellini, L. & Haywood, M. 2010, ApJ, 714, 1096 
Steinmetz, M., Zwitter, T., Siebert, A., et al. 2006, AJ, 132. 1645 
van Leeuwen, F. 2007, A&A, 474, 653 

West, A. A., Morgan, D. P., Bochanski, J. J., et al. 2011, AJ, 141, 
97 

Woolf, V. M., Lepine, S. & Wallerstein, G. 2009, PASP, 121, 117 
Wielen, R. 1977, A&A, 60, 263 

Yanny, B., Rockosi, C, Newberg, H. J., et al. 2009, AJ, 137, 4377 
York, D. G., Adelman, J., Anderson, J. E., et al. 2000, AJ, 120, 
1579 

Zwitter, T., Siebert, A., Munari, U., et al. 2008, AJ, 136, 421 
Zwitter, T., Matijevifi, G., Breddels, M. A., et al. 2010, A&A, 522, 
A54 



